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Abstract 

Intrinsic quantum correlations of weak coherent states are observed between two parties through 
a novel detection scheme, which can be used as a supplement to the existence decoy-state BB84 
and differential phase-shift quantum key distribution (DPS-QKD) protocols. In a proof-of-principle 
experiment, we generate bi-partite correlations of weak coherent states using weak local oscillator 
fields in two spatially separated balanced homodyne detections. We employ nonlinearity of post- 
measurement method to obtain the bi-partite correlations from two single-field interferences at 
individual homodyne measurement. This scheme is then used to demonstrate bits correlations 
between two parties over a distance of 10 km through a transmission fiber. We believe that the 
scheme can add another physical layer of security to these protocols for quantum key distribution. 
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Quantum entanglement and superposition provide secure communication between two 
parties for key generation and information processing. However, entanglement based key 
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generation such as Ekert's protocol [1] is hard to implement in real- world optical fiber net- 
work because bi-partite correlations of entangled photon-pairs are sensitive to loss. There 
has been much interest in quantum key generation of using weak coherent states or highly 
attenuated lasers. Quantum key distribution using weak coherent states, such as coherent 
state differential phase-shift quantum key distribution (DPS-QKD) {2, 3] and decoy-states 
BB84 protocols |4j-l9|, have been proven to be unconditional secure against photon-number 
splitting attack (PNS). The DPS-QKD uses intrinsic correlations between the relative phase 
shifts {0, 7r} of two consecutive pulses to achieve unconditional security between two par- 
ties by constructing equivalent states for the entanglement-based protocol [3j. The decoy 
state quantum key distribution uses intrinsic correlations between the relative mean photon 
numbers of two set of weak coherent states to detect PNS attack in BB84 protocol 4j]. Mean- 
while, Y00 protocol [10] uses intrinsic correlations between phase and mean photon number 
fluctuations of weak coherent states to provide cryptographic service of data encryption 
between two parties. 

Intrinsic quantum correlations of coherent states can be prepared, measured and shared 
between two parties for quantum cryptography. We are motivated to propose a scheme 
based on weak coherent states for generating intrinsic bi-partite correlations as a supplement 
resource to the existence protocols such as coherent state DPS-QKD and decoy state BB84. 

Weak local oscillator (LO) field in a coherent state has been successfully used to directly 
measure bi-partite correlation functions of a two-photon source and violate Bell's inequalities 
using homodyne detection with photon counting 11] . In this work, we employ a weak local 
oscillator field in a coherent state to extract intrinsic correlations of weak coherent states 
between two parties using balanced homodyne measurement. Briefly, we first prepare a weak 
coherent state using a highly attenuated laser at telecom wavelength. The coherent state is 
split by a 50/50 beam splitter and sent to Alice and Bob. Alice and Bob, each has a balanced 
homodyne detection scheme for measuring his/her coherent state with a weak local oscillator 
field. We employ nonlinearity of post-measurement method, i.e., multiply two single- field 
interferences from individual balanced homodyne measurement. Then, the mean- value of 
the multiplied signal provides raw data of field-field correlations of weak coherent states. 
We normalize the raw data with the mean photon numbers of weak coherent state and LO 
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FIG. 1: Experiment setup for observing coherent state bi-partite correlation of weak coherent 
states. 

field to obtain the coherent state bi-partite correlation function (CSBC). It is named so as 
to avoid confusion with bi-partite correlations provided by four Bells states of polarization- 
entangled photon pair. Four types of correlation functions ±cos2(#i ± 62) can be prepared 
by using linear optics devices in Alice or Bob alone, where 9\ and 82 are the projection angles 
of the analyzers at Alice and Bob. This means that Alice can keep her copy of the coherent 
state and send another copy to Bob. By locally changing the relative phases between her 
coherent state and weak local oscillator field, her acts will change the correlation functions 
shared with Bob. Once we establish one of the four correlation functions between Alice and 
Bob over a distance of 10 km through a transmission fiber, we change the phases of the weak 
local oscillator field {0, 7r} for implementing bits correlations between them. 

The weak coherent state |a) and weak LO field in a coherent state \/3) are treated as a 



product state of two independent coherent states \a, (3) in the input of the beam splitter [12 ]. 
The density matrix of the output state from the beam splitter is depended on the integra- 
tion over the phase space of P- representation for the input product state. Since the \ot) 
and are intrinsically correlated from the same laser through the LO phases {0, 2ir}, then 
integration over the phase-spaces of a and (3 can produce the output state that is intrinsi- 



cally entangled [13j | . This is accomplished in our experiment by conducting the mean- value 
measurement of the multiplied signals of the output state. The two mode coherent states 
at the output of the beam splitter, \a + (3)i\a — (3) 2 , can be manipulated by linear phase 
shifters to project out coherent interferences parts \ai(3 2 ) ± \ft1CK2) or laia^) =t \P1P2) plus 
phase-space noises of |a) and \f3). 

A proof-of-principle experiment for the brief description above is shown in Fig.l. A 
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highly attenuated laser at telecom-band wavelength of 1534 nm is used to provide weak 
coherent states and weak local oscillator field. We use a 50/50 beam splitter (BS) to mix 
the horizontally polarized coherent state \a) and the vertically polarized local oscillator field 
The input field operators a and b at the beam splitter are the annihilation operators 
for the coherent state |a) and the LO field \/3), respectively. Using unitary transformation 
matrix of a 50/50 beam splitter, the output modes at the beam splitter are cj = -^(dx) + i6y 
in beam 1 and d 2 = -^(idx) + by in beam 2. A quarter wave plate (QWP) at 45° is 
inserted at beam 1 and beam 2 to transform the linearly polarized states to circularly 
polarized states. Then, after the matrix transformation of a quarter wave plate, the field 
operators 6\ — > c± = + b)x + i(b — a)y) and d\ — >■ di = -^{(ia — b)x + (b + a)y). 

A half- wave plate HWP1 (HWP2) is inserted in beam 1 (2) before a cube polarization 
beam splitter PBS1 (PBS2) to project out the polarization state 6 1 (6 2 ) with unit vector 
e i{2) = cos#i( 2 )X + sm6*i(2)y, respectively. The field operators c\ = (c/ • ei)ei after the 
PBS1 at beam 1 and d 2 = \d 2 ■ e 2 )e 2 after the PBS2 at beam 2 provide the photon number 
operators of the transmitted component of the combined fields as, 

// 1 

C* ci' = - [a*d + b*b + d*be l2dl + b*de~ l2dl ] (1) 

d*"d 2 " = -[ a *a + b*b - d*be i2d2 - b*de" ad2 }. (2) 

Now, using the input weak coherent state \oc) = \ae lr ^ a ) and the input weak LO field in 
a coherent state \0) = \fie l ^ LO ) into Eq. [1] and Eq. [2J the detectors A and B measure the 
transmitted components of the beat intensities, 

(01, 0lo) -> (a, f3\6*"c"\a, /3) = VA [\a\ 2 + \(3\ 2 + \a\\f3\Cos(26 1 + LO - <f> a )] (3) 

I 2 || (62, <!>lo) -> (a, 0\d*"d"\a, 0) = ^[|«| 2 + \/3\ 2 - \a\\(3\Cos{26 2 + LO - <f> tt )\ (4) 

where \a,f3) is the input state of the beam splitter and T)a(vb) is the conversion efficiency 
(Watt — > Current) for detection electronics of A (B). The first two terms are intensities of 
the two coherent states and the last term is the interference term consists of polarization 
angle 6\ (9 2 ), the phases of LO (4>lo) and the weak coherent state (</><*) . The beat intensity 
h\\(9i, (pLo) is anti-correlated to I 2 \\(0 2 , <p LO ) because of the 7r-phase shift induced by the 
50/50 beam splitter. The beat intensity for the reflected signal at the PBS1 (PBS2)is 
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h±(Oi±,<f>Lo) (h±(Q2±,<f>Lo)), where the 9 1± = 9 X + vr/2 (9 2± = 9 2 + 7r/2), respectively. 
Then, the detectors A and B measure the balanced homodyne beat intensities, that are, 

Ai ->■ h\\ - h± = 277A|a||/3| cos(2^ + LO - <t>a) (5) 

/3 2 -> J 2 || - I 2± = -2 VB \a\\(3\ cos(2# 2 + <f) LO - <j> a ) (6) 

plus shot noise at each detection of A\ and B 2 . The individual intensities of the two coherent 
states are subtracted. One can see that the information 2# 1 +0 LO and 29 2 +4> LO are protected 
by quantum phase noise (f> a with phase fluctuation of A(p a > where low mean photon 
number fluctuation associated with phase fluctuation is provided by the weak coherent state. 
The balanced homodyne beat intensities in detectors A and B are then multiplied to obtain, 

A X E 2 = -2r ]A r ]B \a\ 2 \f3\ 2 [cos2(9 1 - 9 2 ) + cos 2{9 1 + 9 2 + <p LO - (7) 

The multiplied balanced-homodyne beat intensities did not provide bi-partite correlation 
function directly. This is predicted because so far we have performed the multiplication 
of single-field interferences obtained from individual balanced homodyne detection at the 
detectors A and B. However, by taking the mean value of this multiplied beat intensities, the 
last term is averaged to zero due to slowly varying local oscillator phase 4>lo from {0, 2tt} 
protected by randomness of quantum phase noise <f> a . Note that A(p a cannot randomly 
provide phase shift — > 2n. And hence, we obtain the expectation value of two detectors or 
coherent state bi-partite correlation function as given by, 

AJEh ->■ (AiB 2 ) = -2r ]A r ]B \a\ 2 \/3\ 2 cos2(9 1 -9 2 ). (8) 

The coherent state bi-partite correlation is protected by the term cos2(#i + 9 2 + 4>lo — (fia), 
which is averaged to zero. In real practice, the product of mean photon numbers |a| 2 |/9| 2 
is obtained by setting the correlation function to its maximum obtainable value, that is, 
Q\ — 9 2 . The raw data of the multiplied beat intensities is then normalized with the product 
of 277^4770 1 a | 2 1/3 1 2 to obtain correlation function — cos2(#i — 9 2 ). As for the use of quantum 
communication between two distant observers, we have to establish the bi-partite correlation 
function and then to implement bits correlations at detectors A and B through the LO 
phases {0, n} without using the post-measurement method. To verify the analysis above, we 
perform systematic studies of the experiment. The balanced homodyne detectors are made 
of two PIN photodiodes (EXT500). We use a piezoelectric transducer (PZT) to ramp the 
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FIG. 2: (a) Shot noises of the weak LO field (0 dBm) and coherent state (-15 dBm), and the 
corresponding beat signals for 9\ = 62 at detector A (b) and detector B (c), and their multiplied 
signal (d). (e) shot noises for the weak LO field (-30 dBm) and the corresponding beat signals for 
9\ = 62 at detector A (f) and detector B (g), and their multiplied beat signal (h). (i) the multiplied 
beat signal when the relative angle 0\ — O2 = ir/2. (square dot) Electronic noise. 

phase of the weak LO beam. We first perform the experiment by using a strong LO field 
and a coherent state with average power of dBm and -15 dBm, respectively. Figure 2 (a) 
shows the spectrum of the shot noise levels of the strong LO field and coherent state, and the 
electronic noise of our detection system. We set the relative angle between Alice's and Bob's 
analyzers as Q\ — 9 2 = 0. The beat intensities at detectors A and B are shown in Fig. 2(b) 
and (c)(inset of Fig. 2(a)). The discontinuity of the beats are due to the ramping of the PZT. 
With these large mean photon fluxes, the interference signals are stable as predicted by the 
coherent states with large mean photon number. The product of the beat intensities is shown 
in Fig.2(d), which is a clear indication of anti-correlation between single- field interferences 
at detectors A and B. Then, we attenuate the laser light to obtain weak LO field and weak 
coherent state with average power of -30 dBm and -30 dBm, respectively. All the average 
optical powers reported in this work are measured just before the PBS1(PBS2). Figure 
2(e) shows the shot noise of weak LO field almost falls on the electronic noise spectrum. 
We observe the beat intensities at detectors A and B as shown in Fig. 2(f) and (g) (inset 
of Fig. 2 (a)) with the interference signals hidden or protected by the shot noises of the 
LO field, quantum phase noise <p a of the weak coherent state due to low mean photon 
number fluctuation, and electronic noises. These are predicted by the Eq. [5] and Eq. [6j In 
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the experiment, the beat intensities at detectors A and B are stored in their computers. 
Then, these raw data are multiplied together as shown in Fig. 2(h). The multiplied beat 
intensity consists of two parts; coherent and noise interferences. The coherent interference 
part contains the term — 2^^|«| 2 |/3| 2 cos 2(9\ — 8 2 ). The noise interference part contains 
the term — 2^775 1 a | 2 1 /3 1 2 cos2(6'i + 9 2 + 4>lo ~ 0a) which is averaged to zero because of the 
periodic of the LO phase, (fiio = {0,27r}, protected by the quantum phase noise (fi a . The 
contribution of short noise is averaged to zero. Since the electronic noise is not completely 
random and present in our measurement method, the noise will create statistical errors in 
the mean-value measurement of the coherent part of the multiplied signal. Note that the 
measurement method is also applied for large mean photon number coherent states as shown 

n 

in Fig.2(d) and also for the mixture of stable and noise fields [14J ]. The aim of the post- 
measurement method is to make sure that weak coherent state for quantum key distribution 
can provide coherent state bi-partite correlation, which can be used as a supplement for the 
decoy state BB84 and coherent state DPS-QKD. From the Fig.2(h), we obtain the product of 
the mean photon number fluxes for the |a| 2 |/3| 2 , where the — cos 2 (#1 — 9 2 ) = — 1 is maximum 
obtainable value for Q\ = 9 2 . Fig. 2(i) shows that the multiplied signal is proportional to 
— cos 2(6^ — 2 ) — 1 when the relative angle is set to Q\ — 8 2 = n/2. We are able to prepare 
four types of bi-partite correlations such as ±cos2(6*i ± 9 2 ) shared between two parties by 
using liner phase shifters on either beam 1 or beam 2. For practical quantum communication, 
Alice can keep the beam 1 and linear phase shifters as highlighted in the box in Fig.l, and 
send out the beam 2 to Bob. Since Alice can change the phases of beam 1 locally, her 
acts will change the coherent state bi-partite correlation function shared with Bob. Fig.3(a) 
shows that the normalized coherent state bi-partite correlation function — cos 2(#i — 6 2 ). For 
each data point, we take 10 shots of the multiplied signal and obtain the average mean-value. 
The error bar is mainly due to the electronic noises. The offset of the relative angle due 
to imperfection of quarter wave-plates has been corrected. For preparing the correlation 
function of cos2(#i + 8 2 ) in Fig. 3(b), we set the A/4 wave-plate in beam 1 to -45° so that 
the beat intensity A\ of Eq. [5]becomes A\ oc — cos(2#i — (0lo — 4>a))- As for the correlation 
function of cos2(#! — 9 2 ) in Fig. 3(c), we insert a A/2 plate at 0° in beam 1 so that the 
minus sign of beat intensity A\ of Eq. [5] is changed to positive sign. Similarly, with the 
A/2 wave-plate at 0° and the A/4 wave-plate at -45° in beam 1, the beat signal A\ of Eq. [5] 
becomes cos(26' 1 — (4>lo ~ 4>a))- Thus, the correlation function of — cos 2(6*i + 6 2 ) is obtained 
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FIG. 3: The coherent state bi-partite correlations for (a) (box) — cos 2(8% — 8 2 ), 

(b) (diamond)cos2(6>i +0 2 ), (c)(cross)cos 2(6»i -82), and (d) (circle) -cos2(0i + 8 2 ). 

as shown in Fig. 3(d). 

After we establish or choose one of coherent state bi-partite correlation functions between 
Alice and Bob, we implement bits correlations between them. To perform this measurement 
for the established correlation function of — cos 2(6^ — 8 2 ), we ramp the PZT to obtain 
one period of interference signal. We reduce the average power of the weak LO field to 
-39 dBm and the average power of the weak coherent state to -39 dBm. The output of 
the balanced homodyne beat intensity at the detector A is directly connected to a lock-in- 
amplifier, where the reference frequency at 168 Hz is obtained from a function generator that 
drives the PZT. We measure quadrature phases of weak coherent state with the step size of 
nit /2 (n=integer) as shown in Fig.4(a) (solid line). Using the same lock-in reference phase in 
the lock-in-amplifier, we measure the quadrature phases of weak coherent state at detector 
B as shown in Fig.4(a)(dashed-line). We have observed the bits correlations between two 
parties for the shared correlation function of — cos2(#! — 8 2 ) as shown in Fig. 4 (a), where 
the positive (negative) quadrature signal is encoded as keys/bits '1' ('0'), respectively. By 
using the same lock-in reference phases, we observe bits correlations for other three types of 
correlation functions — cos2(#! + 8 2 ), cos2(6' 1 + 8 2 ), and cos 2(#i — 8 2 ) as shown in Fig. 4(b), 

(c) and (d), respectively. In real practice, we can establish one of the CSBC for calibrating 
the lock-in reference phases at Alice and Bob. 

We test the scheme by performing bits correlations between two parties over a distance 
of 10 km through a transmission fiber. We couple the beam 2 into the transmission fiber. 
A quarter wave plate and a half wave plate (not shown in Fig.l) are used at the output 
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FIG. 4: Bits correlations for CSBC of ((a)and (e)) - cos 2(0i-0 2 );((b)and(f))- cos 2(0i+0 2 );((c)and 
(g)) cos 2(0i + 2 );((d) and (h))cos 2(0i - 02) , where (e), (f), (g) and (h) are with the 10 km fiber. 
The size of the data point is the error bar for the measurement. 

of the transmission fiber to compensate the birefringence. Since there are losses in the 
couplings and the fiber, we use weak coherent state with average power of -33 dBm and 
weak local oscillator field with average power of -33 dBm before the homodyne detections. 
We have established four types of correlation functions and performed bits correlations for 
each shared correlation function between two parties as shown in Fig. 4(e), (f), (g) and (h). 

As a supplement to the DPF-QKG, the phase of the weak LO field can be randomly 
modulated as {0, tt} at certain frequency. Then, the bits/keys correlations can be realized 
based on the established CSBC shared by both parties. Since the established CSBC is 
normalized with the product of mean photon numbers |a| 2 |/3| 2 , photon number splitting 
attack can be detected by adding a weak LO beam in the decoy state BB84 protocol to 
check the CSBC shared between two parties. The security analysis of the scheme is out of 
the scope of this paper. 

Intrinsic correlations of coherent light field have been utilized to implement entangle 



ment 



15] , Grover search algorithm 
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171 ]. quantum lithography [18| . factoring number 19] 
I I 

and quantum walk [20] through different well-designed interference measurement methods. 
Intrinsic correlations of coherent states do not exhibit non-locality as two-photon source. 
The realization of intrinsic quantum correlation of weak coherent state by using the mea- 
surement method is a first step toward linear-optics quantum computing with weak light 
fields and single-photon source. 
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